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Role of friction-induced torque in stick-slip motion 
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PACS 46 . 55 . +d - Tribology and mechanical contacts 
PACS 46.50.+a - Fracture mechanics, fatigue and cracks 
PACS 81 .40.Pq - Friction, lubrication, and wear 

Abstract. - We present a minimal quasistatic ID model describing the kinematics of the transition 
from static friction to stick-slip motion of a linear elastic block on a rigid plane. We show how 
the kinematics of both the precursors to frictional sliding and the periodic stick-slip motion are 
controlled by the amount of friction-induced torque at the interface. Our model provides a general 
framework to understand and relate a series of recent experimental observations, in particular the 
nucleation location of micro-slip instabilities and the build up of an asymmetric field of real contact 
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Interfacial friction [IH1] plays a major role in seismol- 
O ogy [5] , biology PPS] and nanomechanics |3] . The frictional 
r/3 behavior of a contact interface is controlled by the shear 
k-*^Btrcngth field cr c (x), with x the position in the interface. 
*^As the shear force is increased, a slip region nucleates at 
L_Jhe first point where the shear stress reaches the shear 
strength and grows through the propagation of a micro- 
^Jslip front. Macroscopic sliding starts only after the entire 
^ interface has slipped. This general picture has provided 
OQ^he basis for friction models for decades [7HI5]. Recently, 
("SJ transitions from static to kinetic friction received renewed 
CNl interest due to experiments that allow the local dynam- 
v^ics of frictional interfaces to be directly measured (see e.g. 
<^> |1 4f[T7] ) . Depending on the contact configuration, differ- 
ent kinds of transitions are observed. 

. For contacts between bodies having different shapes (e.g. 
J^j a sphere on a plane) the transition is smooth. As the shear 
force is increased, micro-slip occurs immediately at the pe- 
C3 riphery of the contact where the pressure vanishes, and the 
slip zone quasistatically invades the higher pressure cen- 
tral region [TSl[r7]. For multicontacts (contacts between 
rough solids), this behavior was predicted decades ago by 
Cattaneo and Mindlin using Amontons' law of friction 
(cr c (x)=^p(x) where p is the pressure and fi is the fric- 
tion coefficient), as follows. The distribution of pressure 
p(x) and shear stress er(x) for a non-slipping interface is 
first calculated. In an annular region at the contact's pe- 
riphery, cr(x) > tr c (x), showing that slip must take place. 
Such knowledge of the geometry of the slip region allows 
the stress field for partial slip conditions to be calculated 



|18| . In the following, we will apply a similar procedure to 
planar contacts. 

For planar contacts, er c (x) is expected to be essentially 
homogeneous and micro-slip nucleation should occur at 
random locations due to unavoidable heterogeneities. Sur- 
prisingly, in most experiments, and whatever the way the 
normal and shear loads are applied, micro-slip starts at 
the trailing edge of the contact and propagates dynami- 
cally towards the leading edge [PollTM2~i] . Recent experi- 
ments on multicontact interfaces also show that the onset 
of sliding is preceded by a series of precursors of increasing 
length, which arrest before reaching the interface's leading 
edge [2"Tl[251IM] . These precursors have also been observed 
indirectly in microstructured contacts [22] . 

Two aspects of these findings triggered an active de- 
bate. One concerns the dynamics of the three types of 
micro-slip fronts observed in [T5J[2TJ[33] . This variety has 
recently been reproduced in a ID spring-block model in 
which the interface obeys a complex dynamics described 
by arrays of springs having a distribution of detachment 
force thresholds and a constant delay time for reattache- 
ment [13]. This first aspect will not be further addressed 
here. 

The second aspect, which is the main focus of this Let- 
ter, concerns the kinematics of the precursors, i.e. how 
their triggering force, number and length arc selected. 
Precursors are always found to nucleate near the trailing 
edge. They are accompanied by the growth of an asym- 
metric field of real contact area, with a minimum near 
the trailing edge, which is retained during macroscopic 
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Fig. 1: Sketch of the system geometry (see text for details) for 
X — (dashed line) and X/0 (solid line). 



motion (2TJ[23]- Conversely, changing the way the shear 
[2lJ[23] or normal [24] load is applied modifies the series 
of precursors. To date, no general description includes 
all these kinematic features. However, these observations 
strongly indicate that precursors are highly dependent on 
how both the pressure and shear stress are distributed at 
the interface. 

Experiments on planar contacts can be sorted into two 
groups according to the choice of macroscopic loading: 
the slider is either pushed from its trailing edge side 
[T5J HU [531 [21] or driven from its top by a rigid body 
19,20,22 . In the first configuration, it has been sug- 
gested that slip must nucleate at the trailing edge because 
pushing dramatically increases the shear stress in its vicin- 
ity [T3H23]- I n the second configuration, the shear stress 
is expected to be distributed homogeneously over the in- 
terface, and the reason why nucleation also occurs at the 
trailing edge is still unclear. In the following, we will focus 
on this second, top-driven, configuration because the de- 
scription of the stresses is much simpler than for the first 
(which would in particular include the contact stress field 
around the pushing point). 

In this Letter, we present a minimal ID quasistatic 
model for a sheared planar frictional interface. It involves 
a stress analysis which is inspired by that of Cattaneo and 
Mindlin. By accounting for the first time explicitly for 
the torque that arises whenever the tangential force is not 
applied exactly in the plane of the interface, we reproduce 
most of the above-mentioned kinematic observations. This 
provides a comprehensive picture of the kinematics of the 
transition from static friction to periodic stick-slip motion, 
including precursors. 

Model. We consider the setup sketched in fig. [TJ A fric- 
tional interface is formed by pressing a linear elastic block 
(Young's modulus E, thickness t, width w and length L, 
perfectly bonded to a rigid plate) on a horizontal plane. 
The normal load N, applied symmetrically with respect 
to the contact center, induces no tilt of the plate. The 
tangential displacement of the plate, X , is prescribed at a 
height H with respect to the frictional interface, giving rise 
to a friction force F. We assume Amontons' rigid-plastic 
law of friction, with a static friction coefficient fx s and a 
kinetic friction coefficient ^ < Ms- We therefore neglect 
any effect of the tangential stiffness of the multicontact 
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Fig. 2: Calculated stress fields p(x) (open symbols), and cr(x) 
(full symbols) for F=0 (circles), 0.41JV (squares) and 0.827V 
(diamonds). Stresses are made non-dimensional by dividing 
them by the average pressure. For clarity, only 1 data point 
out of 9 is plotted. These fields are compared to the fields 
assumed in the model (dashed lines). 



interface, which would induce deviations to Amontons in 
contacts where stick and slip regions coexist [25] . We also 
neglect state and rate effects [3] . We place the origin of x 
at the center of the contact, which extends between 
The problem is made dimensionlcss by expressing coordi- 
nates in units of L, forces in units of [i s N and stresses in 
units of ^j^- Dimensionless quantities bear a tilde. All 
physical quantities can be expressed in terms of only two 
dimensionless control parameters: r = *jf- and g = 

In order to guide our assumptions about the stress dis- 
tributions we first performed simple, plane strain, finite 
element calculations (software Castem 2007) on the con- 
figuration sketched in fig. [TJ The boundary conditions 
are the following: no displacement is allowed at the con- 
tact interface ; the linear elastic block (L=40cm, i=2mm, 
i?=2.10 3 Pa, Poisson's ratio 0.4) is perfectly bonded to the 
rigid plate (75=2.10 15 Pa, Poisson's ratio 0.4) ; the normal 
displacement of the rigid plate is prescribed at the center 
of its top surface ; the tangential displacement of the rigid 
plate is prescribed on a point of its left side at a heigth 
H=4cm. We used a regular mesh size of 1mm and QUA4 
elements. Figure [2] shows the pressure and shear stress 
fields over the contact interface for normal loading only 
and for two different shear forces applied. For F—0, the 
shear stress is zero and the pressure field is a constant. For 
F t^O, the shear stress increases homogeneously over the 
whole contact, whereas the pressure develops an asymme- 
try which, to a very good approximation, is linear with 
x. We find deviations at both contact extremities due to 
edge effects, which are significant over distances ~ t for 
the pressure and ~ 3t for the shear stress. Based on these 
preliminary calculations, we develop the following model. 
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Fig. 3: Stress fields fi s p(x) (solid line), fidi>(S:) (dashed line) and a(x) (dotted line), (a) F = 0: all fields are homogeneous and 
cr(x) = 0. (b) F = F±: fj, s p = a = F\ at x = —1/2 (black disk). Propagation stops at xi (open disk) where fidft — 5" = F°. 
Vertical arrows show the stress relaxation on the segment [— 1/2; xi]. The grey surface represents the relaxed force F*. (c) 
Same as (b) but for F = F2. 



We neglect any edge effect like e.g. a divergence of the 
stresses at the border of the punch-like contact Q2] , an as- 
sumption that is increasingly valid as t/L decreases. These 
edge effects would be symmetric with respect to x — and 
hence could not account alone for the observed asymme- 
try. Under these conditions, for F = the pressure field 
p(x) = J- is therefore constant (fig. [3Ja)). Moreover, be- 
tween two micro-slip events, the tangential displacement 
of the top surface of the block is homogeneous, so that the 
shear stress field &{x) increases homogeneously (fig.[3tb)). 

For F ^ 0, a torque FH is applied at the interface, 
which can only be balanced by an asymmetry of the pres- 
sure field. The elastic block is confined between two rigid 
planes, one of which (the plate) undergoes a slight tilt. 
This yields a linear spatial distribution for the normal 
compression, which is assumed to translate into a linear 
pressure field p(x, F) = — + j^Fx, which agrees very well 
with the one calculated above (see fig. [2]). This assump- 
tion is increasingly valid with decreasing t/L, and would 
increasingly break for an increasing compliance of the driv- 
ing plate. We assume that this linear form remains true 
even after micro-slips, i.e. the pressure field, just like the 
torque, depends only on the total friction force F and not 
on the distribution of shear stress over the contact. This 
is similar to the classical Goodman assumption in contact 
mechanics which is here increasingly valid with 

decreasing t/L. 

At a certain force F = F\ the local static slip threshold 
cr(x) = fi s p(x, F) is reached for the first time at the trailing 
edge x = — h (if g > 0), where the pressure is minimum 
(black disk in fig. E^b)), and we find that: 

Fi = T-^- (1) 

1 + g/r 

A micro-slip front nucleates at the trailing edge, turning 
the problem into the one of an interfacial shear crack. 



This crack is unstable in the sense of the Griffith ener- 
getic criterion (sec Appendix) over virtually all x such that 
d'(x) £ \pLdP; HsP\- The micro-slip front therefore prop- 
agates towards the leading edge. By assuming that the 
friction force relaxation associated with the micro-slip oc- 
curs only after propagation is over, we find that the front 
arrests at point x\ such that a{x\) = fj,dp(xi, i*i) (open 
disk in fig.ES»). 

In all the slipped region, the shear stress drops to 
fj,dp(x,Fi) (arrows in fig. [3jb)) so that the shear stress 
field is now a(x) = r + 2gF\x for x € [— 1/2; ii] and 
a(x) = F\ elsewhere. We therefore neglect the exten- 
sion ~ t over which the slope discontinuity of the shear 
stress at iii is regularized. The relaxed friction force F* is 
the integral of this relaxed shear stress field (grey area in 
fig-ElJb)). The subsequent evolution for increasing F is a 
series of micro-slip events following the same scenario (see 
fig. G2c) for a sketch of the second event). 

We now derive the iteration formulae for the successive 
values of Fi for such precursors, the corresponding arrest 
positions Xi and the relaxed forces F*. We introduce the 
elastic force F c = KX, where K is the effective stiffness 
of the elastic block. F c is the tangential force that would 
have been required to move the rigid plate in the absence of 
partial relaxations related to the micro-slip events. Before 
macroscopic motion, F c reduces to a(x = 1/2) (see fig. [3]). 
The value Ff at the onset of the i-th precursor is: 

i-i 

F * = F + E - F j) ( p i = ( 2 ) 
3=1 

The arrest point for the i-th event is then given by 
fi d p{xi,Fi) = F? (see fig. [1(b) and (c)), yielding: 

Xl = (P? - r) /2F l9 . (3) 
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Knowing Fi, F? and Xi, the relaxed force F* is derived 
by integrating the shear stress field just after the precur- 
sor. One shows graphically (fig. [3l grey areas) that: 



F* = F c 



1 

- {Xi 

2 



(ii + 1/2) Ff-rll-gFi/r 



(4) 



We finally calculate the force at the next event. The static 
threshold is always reached first at x = —i. Then F + i 
verifies /j, d p{-^, Fi)+F i+1 -F* = /Lt s p(-|, F i+1 ), yielding: 



F, 



i+i 



Ft+gFi 



1 + g/r 



(5) 



The evolution of the system can be solved itcratively with 
Eq. Q]as the starting point. At each step, F* , F i+ i, Ff +l 
and Xi+\ are computed successively. This procedure has 
to be slightly modified as soon as x n +\ > 1/2, i.e. when 
the (n + l)-th event first reaches the leading edge after 
n precursors. At this point, the shear stress relaxes to 
a{x) = (j,dp(x) everywhere, i.e. F* +1 = r. The subsequent 
evolution is then obtained by simply replacing Eq. [3] by 
Xi = 1/2 and Eq. 0]by F* — r in the iteration. 

Results. Figure [?] shows how the friction force F typ- 
ically evolves with the elastic force F c (which is propor- 
tional to the prescribed displacement X). This loading 
curve exhibits an initial elastic regime at low forces. It 
is followed by a series of n relaxations corresponding to n 
precursors of increasing length. The system then enters a 
macroscopic stick-slip motion regime (each slip event in- 
volves the whole interface), in which the relaxed force is 
always r. Stick-slip progressively becomes periodic, reach- 
ing a maximum friction force of F max = 1+ ±^ 1 _ r ^ (calcu- 
lated by imposing F i+ i = Fi in Eq. [5]). 

In practice, the two control parameters can only be var- 
ied in the range 0<r<l(r>lis not physical because 



fi s > fid) and < g < 1 (p > everywhere in the con- 
tact in the absence of adhesive forces, which translates 
into the condition that, for all F < F maX7 p(— 1/2, F) > 0, 
yielding r — gF max > and eventually g < 1). Figure [5] 
shows the number of precursors n over the whole accessible 
parameter space, n increases with both r and g. In par- 
ticular, the larger the torque, the more precursors. Below 
the line g = (1 — r)/2, the extension of the first precursor 
X\ is longer than the contact size, so that no precursor 
is observed. In particular, no precursor can occur at a 
torque-free interface (H = 0, i.e. g = 0). 

Discussion. We built a minimal model accounting for 
the rich kinematics of the transition from static to kinetic 
friction. The key ingredient is the increasing asymmetry 
induced by the growing tangential force, as soon as the lat- 
ter is not applied exactly in the plane of the contact inter- 
face. We emphasize that this effect is different in essence 
from any time-invariant pressure asymmetry that would 
be due to an asymmetric normal loading, like in [24] . Our 
strongest assumption is that the micro-slip front is so fast 
that the associated force relaxation occurs only after its 
arrest. This is a hidden assumption about the propaga- 
tion dynamics, which is otherwise explicitly beyond the 
scope of this quasistatic model. In this respect, our model 
is a crude description of a real frictional interface. Yet, 
we believe that it captures important generic aspects of 
sheared frictional interfaces. 

An important component of the present model is the 
build up of a friction-induced torque. We emphasize that 
in typical situations, a significant torque corresponds only 
to a very slight tilt angle. For instance, for small defor- 
mations F/wLE < 10~ 3 , a force applied at the bottom 
of the rigid plate H — t and an aspect ratio t/L < 1CP 1 
yield a < 1.2 • 10~ 4 rad. The asymmetry of the contact is 
therefore very difficult to avoid experimentally. 

Even for an initially homogeneous interface, tangential 
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loading produces an increasing contact pressure asymme- 
try. In practice, the tangential force is usually applied 
above the interface (H > 0), making the pressure lower 
at the trailing edge. Under these conditions, any friction 
law that prescribes a monotonic increase in shear strength 
with increasing local pressure implies that the lowest shear 
strength is at the trailing edge. In practice most friction 
laws belong to this category. Here, we considered the clas- 
sical Amontons' law, which is known to be valid for macro- 
scopic multicontact interfaces [IJ[27]. The smooth glass- 
on-gel interfaces used in [TjJJ[5D] , although a very different 
system, also obey a pressure-increasing shear strength law, 
which is a modified version of Amontons' law including ad- 
hesive forces. The pressure dependence of shear strength 
of the microstructured contacts studied in |22j was not 
reported but we speculate that it also follows an increas- 
ing trend. Our model therefore successfully explains why 
micro-slip fronts always occur at the trailing edge in both 
systems [TS1[2D1[22] > which are top-driven. 

In the side-pushed multicontacts studied in [2T1I231I24] . 
we also expect the torque to decrease the pressure at the 
trailing edge. This is confirmed by measurements of the 
distribution of real contact area just after the successive 
precursors and after each slip in the periodic stick-slip 
regime [2TJ[23]. In a multicontact interface, the real con- 
tact area is well-known to be proportional to the applied 
contact pressure [2"TI[2"5] . Our model is therefore consistent 
with their observation of a real contact area asymmetry 
which increases with F, and which is retained and stable 
in the macroscopic stick-slip regime, with a minimum near 
the trailing edge. We emphasize that further quantitative 
comparison between our model and the measurements for 
the number and length of the precursors in [2TJ[331[2I] is 
not possible because the two loading systems are very dif- 
ferent. In particular, we believe that the huge shear stress 
increase near the pushing point will likely play the pri- 
mary role in reducing the shear strength at the trailing 
edge, and dominate the torque effect. We also emphasize 
that, even in top-driven systems [THHM!!!!] , quantitative 
comparison would require knowledge of the effective height 
H of the applied friction force, which is seldom provided 
in the literature. We therefore urge authors to provide, in 
the future, all the details about the loading configuration 
that are necessary to assess the level of torque, like in [2T)] 
where H = 0. 

An implication of the fact that micro-slip very gener- 
ically nucleates at the trailing edge of a contact implies 
that the interface will be in a compressed state as soon as 
the first precursor occurs. This strongly suggests that the 
results of Bcnnewitz et al |22] , which measure such com- 
pression, arc valid over a much broader range of systems 
than their microstructured PDMS on glass interface. 

Other results of the model arc of general interest for 
frictional systems. First, they suggest that, from a kine- 
matic point of view, macroscopic stick-slip is not different 
from the preceding precursors, the latter being merely de- 
fined by their limited extent over the interface. Second, 



we show that F max /N is always smaller than i.e. the 
macroscopic friction coefficient is always lower than the 
local friction coefficient, which governs the onset of the 
micro-slip instability. The macroscopic friction coefficient 
is moreover dependent on g, i.e. on the details of the load- 
ing system, suggesting that care has to be taken when a 
friction coefficient is measured using a macroscopic exper- 
iment. 

To conclude, we have shown that in virtually all friction 
experiments, slip will occur at an asymmetrically loaded 
interface due to a friction-induced torque. In particular, 
the pressure is minimum near the trailing edge, where all 
micro-slip instabilities nucleate. In top-driven systems, 
like the one studied here, the kinematics of the transition 
from static to stick-slip friction is then dominated by this 
macroscopic, system-dependent, asymmetry rather than 
by small scale heterogeneities of the shear strength. This 
is in striking contrast with most friction models, which 
assume a homogeneously loaded interface bearing a rela- 
tively small disorder (see e.g. [7H9lll3j'l. In this respect, 
the present work provides a framework for the extension 
of these models to take into account a macroscopic shear 
strength asymmetry. 

Appendix. — Here we show, using an energetic crite- 
rion, that a crack tip at the frictional interface is unstable 
at virtually all points x satisfying the following stress cri- 
terion <j(x) > f±dp{x). We apply the Griffith energetic 
criterion which states that an existing crack will propa- 
gate if G(x) > T c (x), where G(x) is the energy release 
rate at point x and T c (x) is the facture energy at point x. 
To do this, one needs to regularize Amontons' law of fric- 
tion by describing how the shear stress drops from /i s to 
fid- The easiest way, which is classically used in dynamic 
crack simulations (see e.g. [5"I5U] 1 is to assume a linear slip 
weakening: a drops linearly from fi s p to fidP over a critical 
weakening distance D c . Then T c (x) = (^ a ~^ d )- Pc p(x). 

When slip occurs at x, the shear stress is relaxed by the 
amount (n(x) — fid)p(x), where fi(x) = a(x)/p(x). The 
amplitude of the corresponding slip is S(x) = t(fx(x) — 
Hd)^§^. Neglecting the variation of stress during the first 
slip over the (microscopic) length D c , the stress during all 
the slip is fj,dp{x). Therefore G(x) ~ fid(fi{x) — fid)-^p 2 {x). 

The condition G > T r therefore reads 



fx(x) - fl d 



> 



E D cl i s 



fJ-d 



p(x) t 2fi d 



(6) 



D c and t being respectively a microscopic and a macro- 
scopic length scale, D c /t << 1. This means that the crack 
will stop only for very low values of fj,(x) — fid- We there- 
fore consider that a precursor practically propagates until 
the point where o{x) — fidp(x). 
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